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The standard, scale-invariant, inflationary perturbation spectrum will be modified if the effects of 
trans-Planckian physics are incorporated into the dynamics of the matter field in a phenomenological 
manner, say, by the modification of the dispersion relation. The spectrum also changes if we retain 
the standard dynamics but modify the initial quantum state of the matter field. We show that, given 
any spectrum of perturbations, it is possible to choose a class of initial quantum states which can 
exactly reproduce this spectrum with the standard dynamics. We provide an explicit construction 
of the quantum state which will produce the given spectrum. We find that the various modified 
spectra that have been recently obtained from 'trans-Planckian considerations' can be constructed 
from suitable squeezed states above the Bunch-Davies vacuum in the standard theory. Hence, the 
CMB observations can, at most, be useful in determining the initial state of the matter field in the 
standard theory, but it can not help us to discriminate between the various Planck scale models of 
matter fields. We study the problem in the Schrodinger picture, clarify various conceptual issues 
and determine the criterion for negligible back reaction due to modified initial conditions. 

PACS numbers: 98.80.Cq, 04.62.-|-v 



I. THE TANGLED WEB 

The inflationary scenario [1-3] is at present the most attractive paradigm for generating the initial small scale 
inhomogeneities [4, 5]. These perturbations leave their imprints as anisotropics in the Cosmic Microwave Background 
(CMB) [6] and later evolve, due to gravitational instability, into the large-scale structures that we see around us today. 
While there is no natural particle physics candidate for generating the inflationary phase, a single (or a few) 'inflaton' 
field(s), with fine-tuned designer couplings, is (are) often introduced in order to reproduce the observed magnitude 
and shape of the perturbation spectrum. 

In many of the models of inflation [3] , the period of acceleration lasts sufficiently long so that length scales that are 
of cosmological interest today would have emerged from sub-Planckian length scales at the beginning of inflation. This 
suggests that physics at the very high energy scales can, in principle, modify the primordial perturbation spectrum [7] 
and these modifications can — in turn — leave their signatures on the CMB [8] . This has led to a considerable interest in 
understanding the effects of Planck scale physics on the inflationary perturbation spectrum [7]-[20] and the CMB [21- 
23]. 

Metric fluctuations during the inflationary epoch can be modeled by a quantized, massless and minimally coupled 
scalar field [4, 5, 24]. In the absence of a workable quantum theory of gravity, the Planck scale effects on the 
perturbation spectrum have to be studied by phenomenologically modifying the dynamics of the scalar field to take into 
account the quantum gravitational effects (for an early attempt in this direction, see Ref. [7] ) . The high energy models 
of the quantum scalar field that have been popular in the literature either introduce new features in the dispersion 
relation [8-12] or modify the standard uncertainty principle [13, 14] or assume that the spacetime coordinates are 
non-commutative [15]. (For other approaches, see Refs. [16-19] and, for a recent review of many of these approaches, 
see Ref. [20]). Some of these models have been utilized to evaluate not only the Planckian corrections to the standard, 
scale- invariant perturbation spectrum, but also the resulting signatures on the CMB [21-23]. This suggests the 
possibility that sufficiently accurate measurements of the CMB anisotropics can help us understand physics beyond 
the Planck scale. 

There is, however, one serious difficulty with this approach, which we shall now briefly describe. 
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We begin by noticing that, in a Friedmann universe, each mode qi^ of the scalar field, labeled by the wave vector k, 
evolves as an independent oscillator with time dependent parameters that are related to the expansion factor a{t). 
Given the quantum state '^k[9k, ii(k)] for the mode gk at a time ti(k), one can obtain the state at a later time t by 

V'k(g'k,i) = j dq]^K[qk,t; qk,ti{k)]ip]^[q]^,ti{\i)], (1) 

where K{q]f^,t; q)i,ti) is the (path integral) kernel for an oscillator with time dependent parameters, which can be 
written down in terms of the classical solution (see, for e.g., Ref. [25]). To keep our discussion general, we have 
allowed for the possibility that the initial quantum state of each of the oscillators is specified at different times so that 
ti can depend on k. For instance, one may choose to specify the quantum state of each oscillator when the proper 
wavelength of that mode is equal to the Planck length [16, 26]. Of course, when the initial state of all the modes are 
specified at a given time, ti(k) = ti will be independent of k. The dynamics of the system is completely specified by 
the kernel K. 

It is obvious that the mathematics of a free scalar field in a Friedmann universe is trivial and is no more complicated 
than that of an oscillator with time dependent parameters. The perennial interest in this problem (allowing so many 
e-prints to be written!) arises from two conceptual issues: 

1. To make any predictions, we need to know 'ip]i{qk,t), which, in turn, requires knowing ijjk[qk,ti(k)]. We have 
absolutely no idea what to use for tpk[qk,ti{k)] and so different choices (usually called 'vacuum states', which is 
only a manner of speaking) can be investigated. 

2. There are infinite number of such oscillators, leading to the standard, unresolved, issues of regularization in 
quantum field theory. 

To make any progress, we need to make an assumption regarding V'k[3k, UO^)] and our results are only as valid as this 
assumption. 

If we now further modify the dynamics (due to a phenomenological input regarding trans-Planckian physics), we 
will be changing the form of the kernel K. But, since wc can only observe the integrated effect of K and V''k['?k- ii(k)], 
the observations can tell us something about the K (and trans-Planckian physics) only if we assume something about 
"(Aki^k, ii(k)]. The usual assumption is to consider the initial quantum state to be the Bunch-Davies vacuum [27], but 
it is only an assumption. The crucial question is whether the effect of trans-Planckian physics can be mimicked by a 
different choice of the initial state other than the Bunch-Davies vacuum. 

We shall show that any (modified) spectrum of fluctuations can be obtained from a suitably chosen initial state 
V'k[5k,ii(k)], which will prove to be a squeezed state above the Bunch-Davies vacuum in the standard theory. We 
shall provide an explicit construction of the state for any given spectrum of perturbations that is observed. So, if some 
specific deviation from the standard scale invariant spectrum is seen in the CMB, a conservative interpretation will 
be to attribute it to a deviation from the standard initial state of the theory. Unless this possibility is ruled out, one 
cannot claim that the observation supports, say, a particular model of trans-Planckian phenomenology. Motivated by 
this result, we argue that the CMB can at most help us identify the quantum state of the scalar field in the standard 
theory, but it can not aid us in discriminating between the various Planck scale models of matter fields. 

The remainder of this paper is organized as follows. In Section II, we set up the formalism and study the evolution 
of a Gaussian quantum state in a Friedmann universe. We apply this formalism to power law inflation in Section III. 
In Section IV, we show that any modified spectrum can be reproduced from a suitable squeezed state above the 
Bunch-Davies vacuum in the standard theory. We explicitly discuss four modified spectra that have recently been 
considered in the literature. In Section V, we evaluate the energy density in these excited states and examine whether 
these modified spectra also lead to a large back reaction on the inflating background. Finally, in Section VI, we 
conclude with a discussion on the wider implications of our analysis. 

Our conventions and notations are as follows. We shall set h = c = 1 and the metric signature we shall adopt is 
(+,-,-,-)• 

II. EVOLUTION OF THE QUANTUM STATE: GENERAL FORMALISM 

Consider a flat Friedmann universe described by the line-clement 

= _ ^2^^^ ^^2 ^ ^2^^^ ^^^2 _ ^^2^ ^ ^2) 

where t is the cosmic time, a{t) is the scale factor and rj denotes the conformal time with drj = dt/a(t). The scalar as 
well as the tensor perturbations during the inflationary epoch can be modeled by a massless and minimally coupled. 
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real scalar field, say, governed by the action [1, 2, 4, 5, 24, 28] 



= i / dSv^S^^S"*. (3) 
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The homogeneity and isotropy of the Friedmann metric (2) allows us to decompose the scalar field $ as 

(^?Vr/)e^'^- (4) 

The is complex and for each k gives two degrees of freedom in the real and imaginary parts of = At + i-Bk- 
But since $ is a real scalar field, we can relate the variables for k to that for — k by AJ^ = A_k, = — -B_k and only 
half the modes are independent degrees of freedom. Therefore, we can work with new set of the real modes ^k for all 
values of k with suitable redefinition, say, by taking gk = for one half of k vectors and g_k = -Bk for the other 
half. In terms of real variables gk, the action (3) can be expressed as follows: 

S[^ = \ Jd'kJdrja'{qQ-eql), (5) 

where the primes denote differentiation with respect to the conformal time ry and k = |k|. The action (5) describes a 
collection of independent oscillators with time dependent mass and frequencies k. (It is sometimes useful to keep 
track of the real and imaginary parts of the gk separately. In our case, this is unnecessary and we can define our 
system in terms of real gk- Our approach is completely equivalent to the conventional one. Also, note that, in the 
literature, one usually finds the Fourier decomposition in Eq. (4) expressed in terms of (qu/a) rather than with just 
gk- Such a decomposition will lead to oscillators that have a unit mass, but a time- dependent frequency, say, w^, 
which can become negative at super Hubble scales. In our description — which is, again, equivalent to the conventional 
one — the mass varies as a^, but the frequency is constant.) 

In the Schrodinger picture, the scalar field $ can be quantized by quantizing each independent oscillator gk- The 
Hamiltonian corresponding to the k-th oscillator is given by 



Hk=^ + ^a^fc^9^, (6) 



where pk is the momentum conjugate to the coordinate q^. Therefore, each of the oscillators satisfy the Schrodinger 
equation 

and the complete quantum state of the field is described by a wave function that is a product of ipu. for all k. Equiv- 
alently, the time evolution of the wave function can be described by Eq. (1) in terms of the kernel K [<7k, t; q^^, ti(k)]. 

As we do not expect a large scale, spatially inhomogeneous classical scalar field to be present in the universe, it is 
conventional to assume that the expectation value (V'k|^k|V'k) vanishes in the quantum state of the field for k ^ 0. 
(Since (V'kl^klV'k) satisfies the classical equations of motion, this condition can be satisfied at all times if suitable 
initial conditions are imposed at an early epoch.) When the mean value vanishes, the power spectrum as well as the 
statistical properties of the perturbations are completely characterized by the two-point functions of the quantum 
field. Therefore, the power spectrum of the perturbations per logarithmic interval, viz. [A:^P$(A:)], is given by (see, 
for instance, Ref. [28]) 

3 ^ 

k'P^{k) = ^ J dqk\i^k{qk,7i)Wl. (8) 



Though this result is fairly well established in literature, there are a couple of subtleties we would like to mention. 
The quantity on the right hand side of the above equation depends on the time r) and one needs to settle at what 

epoch it has to be evaluated. In classical perturbation theory, one evaluates the perturbation at Hubble exit, i.e. 
when the physical wavelength {k/a)~^ of the mode corresponding to the wavenumber k is comparable to the Hubble 
radius H~^, where H = (a'/a^). In other words, the spectrum is to be evaluated when, say, (k/a) = (zH), where z 
is a number of the order of unity. This is a thumb rule which accounts for the differences in the evolutionary history 
of the mode when its proper wavelength is smaller than the Hubble radius as compared to the situation when it is 
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larger than the Hubble radius. There is no simple way of deciding whether one should evaluate the expression when 
z = 1 or when, say, z = (27r). In the literature, one also finds the perturbation spectrum evaluated at super Hubble 
scales (i.e. when {k/aH) 0) which corresponds to the limit z — > 0. Even for the simplest case of exponential 
inflation with the Bunch-Davies vacuum for the initial state, these results differ by a numerical factor. In this case, 
we will find that, [k^ Pi>(/c)] = [C(z) {W^ /2it'^)], where C{z) = [(1 + z'^)/2] (see Eq. (47) below). It is sometimes 
claimed erroneously in the literature, that evaluating the spectrum as z ^ leads to the same result as evaluating it 
at z = 1. (See, for instance, Ref. [2], pp. 182 183. Notice that, in this reference, Eq. (7.87) is wrong by a factor 2 and 
the claim after Eq. (7.98) that Eqs. (7.98) and (7.87) are "in agreement" is incorrect. This can be trivially verified 
from equations (7.96) and (7.98).) In this particular case (i.e. exponential inflation with the Bmich-Davies vacmim 
for the initial state), this discrepancy is not of great importance since it only changes the amplitude by a numerical 
factor, since C(l) = 1, C(27r) = [(l + 47r^) /2] and C(0) = (1/2). But when one considers the power spectrum in a 
more general context, these choices will lead to a more complicated difference, as we shall see. Of course, if one first 
approximates the wave function i/;k ((/k,??) with an assumption such as z — ^ (when the expressions do simplify) and 
then evaluate it at z ~ 1, one is being inconsistent. This may sound rather elementary but we were surprised to find 
papers in the literature which do this. Any result which crucially depends on one specific choice for z in computing 
the power spectrum is suspect. In what follows, we shall usually assume that the expressions are evaluated for z = 1 
(i.e. when {k/a) = H) when the choice of z is not of much consequence and will comment on the results which depend 
crucially on this choice. 

Let us now consider the problem of determining the wave function i/'k (^k, v)- For a time dependent oscillator, there 
is no concept of a unique ground ('vacuum') state unless the parameters describing the oscillator go to a constant value 
asymptotically — which, in general, it does not, for the Friedmann universe. There is, however, a class of solutions to 
the time dependent harmonic oscillator which are form invariant in the sense that the dependence of the wave 
hmction ?/'k('?k, if) is the same at all 77. One can show that the most general state having this property is an exponential 
of a quadratic function of and such states have been extensively investigated in the past in different contexts (see, 
for instance, Refs. [29, 30]). When (V'kl^klV'k) = 0, the mean value of the Gaussian vanishes and the quantum state 
of the mode can be described by the wave function [28-30] 

V'k (^k, rf) = Nk{ri) exp - [Rk{ri) ql\ , (9) 

where Nk{rj) and Rkiv) are complex quantities. The normalization condition on the wave function then relates Nk 
and Rk as follows: 

\Nkf=(?±±^y'\ (10) 
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Therefore, the only non trivial aspect of the quantum state is encoded in the time dependence of the function 
Rkiv)- It can be shown that (for details, see Appendix A), if we introduce a function Hkiv) through the relation 
Rk = — {ia^/2j (/Ufe//Ufe), then Hk satisfies the differential equation 

^l'l + 2-^^i + enk = o (11) 

a 

which is the same as the classical equation of motion satisfied by the oscillator variable qi^. We find that the pertur- 
bation spectrum (8) corresponding to the wave function (9) can be expressed as (see Appendix A) 

where W{k) is a fc-dependent constant determined by the Wronskian condition for jik [cf. Eq. (A5)]. 

Since the differential equation (11) has real coefficients, if is a solution, so is and the general solution is a 
linear superposition of the form, say, fik — [-^{k) Sk + B{k) si]. The quantum state ■!/'k((?k, v)j however, depends only 
on Rk which is independent of the overall scaling of fik- This feature translates into the power spectrum (12) as well; 
a global scaling of Hk also changes the Wronskian W{k) leaving [|/Xfc|^/W(A:)] invariant. Hence, we can ignore the 
overall scaling in fik- We shall set A{k) to unity and choose the Wronskian W{k) to be 

W{k) = l-\Bik)\'^. (13) 



Then, the power spectrum (12) reduces to 

27r2 I 1- l'B(fc)| 
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(In the Heiscnbcrg picture, one usually introduces the Bogoliubov coefficients a(k) and /?(fc). These coefficients are 
related to B(k) through the relation: B{k) — [/3{k) / a{k)]. In terms of a{k) and /3(fc), the Wronskian condition (13) 
reduces to the standard normalization condition, viz. (|Q!(fc)p — |/3(fc)|^) = 1.) The dynamics imposes no restrictions 
on the form of B{k), which is related to the choice of the quantum state at any specified time. The spectrum evidently 
depends on B{k) and with a suitable choice for B{k) we can generate any, given, reasonable spectrum. 



A. Squeezing and instantaneous particle content of the quantum state 

We shall see concrete examples of this result and its consequences for the case of power law inflation in the following 
sections. But, before we do that, let us try and understand what the wave function (9) implies. The wave function (9), 
in general, describes what is referred to in the literature as a squeezed state (sec, for e.g., Rcfs. [30]). Squeezed states 
for a given mode gk are described by two parameters, say, and ipk and the quantity Rk can be related to these two 
parameters as follows [30]: 

D _(^°^\ ( coshrfc + e^"^" sinhrfc ^ 



2 J \ cosh Tk — 6^"'''= sinh rk 

This, however, does not lead to any deeper insight in this particular case. 

An alternative procedure, which is physically better motivated, is to compare the quantum state (9) with the 
instantaneous ground state at any given time. We recall that the oscillators have the frequency k and a time 
dependent mass a^irj). The ground state of a harmonic oscillator with the frequency k and mass a^{r]) will depend 
on as cxp — [(fc/2) q^]. Suppose we are given a specific quantum state determined by the function Rk{'r]). We 
can attempt to solve the equation Rkiv) — [^^ ^'^ iv) / for rj determining a function, say, fj{k). If such a real solution 
can be obtained, then we can interpret our quantum state as the ground state for the mode qk at the moment of 
time f]{k). In general, however, this will not be possible since Rk{r]) is complex. (One exception is in the adiabatic 
approximation, in which the mode Hk is chosen such that that Rk « (fca^/2) so that this condition is identically 
satisfied; as to be expected, the state evolves as an adiabatic vacuum in this case.) 

More generally, one can expand our quantum state in terms of the instantaneous states of the harmonic oscillator at 
any given time r], thereby identifying its particle content. At any given time r], the instantaneous ground and excited 
states of the oscillators g'k can be described by the wave functions 



0„ (gk,r/) = (2"n!) 



, - . 1/4 

-1/2 



{~^^ ^niykaqkj exp- ^^0^ ql + iEnT]^ , (16) 

where Hn are the Hermite polynomials, En = [{n + (1/2)) k] and n = 0, 1, 2, We can now decompose the actual 

wave function ^k (9k, v) in terms of the above instantaneous wave functions as follows: 

oo 

tpk {qu, ??) = XI ^"f^^' ^" v) e'-^"'', (17) 
where the decomposition amplitude c„ is given by the integral 

oo 

Cn{k,v) = j dqk tpk {qk,v) K (9k,??) • (18) 

oo 

On evaluating this integral, we find that the amplitude c„ corresponding to the odd n's vanish, while the amplitude 

for the even n's are given by 



C,.»,(*..) = A. (^) 



n, (19) 



where 
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is a n-independent normalization and 



(21) 



The probability Pk{n) for our quantum state to be in the (2n)-th excited state of the instantaneous harmonic oscillator 
mode can be thought of as the probability for existence of n pairs of particles at the time rj. This is given by 



Pfc(n)= |c(2„)(fc,r?)|' =Pfe(0) 



(2n) 



nP J \ 4 

and the generating function Qk{o) for this pair creation probability can be expressed in closed form as follows: 



afe(a) = ^Pfc(n)M" 



Pfe(O) 



n=0 



(i-<7|rfe|2)i/2 



1 



i-a|rJ2 



1/2 



(22) 



(23) 



The last equation follows from explicitly calculating Pfe(O) or more simply by noticing that Pk{n) is normalized 
and hence Gk{^) = 1- Given this generating function one can compute various moments of the created particles. In 
particular, the mean number of particles, which are present at the time r] (obtained by doubling the mean number of 
pairs) is given by 



(24) 



These equations exhibit the time dependent particle content of our quantum state and it can be computed once the 
function iikiv) is specified. Though, (n^) can not be interpreted as particles with respect to the in- vacuum at late 
times, it is related in a simple manner to the energy density. The expectation value of the Hamiltonian operator 
corresponding to the oscillator gk, say, £k, can be evaluated using the wave function (9). We obtain that 



£k = 



{Wkf + k'\iik\ 



and, on using the expressions (21), (24) and (25), we find that £k and (rife) are related as follows: 

1^ 



{nk} 



k. 



The energy density of the quantum scalar field is then given by 

oo oc 


We will require these results while discussing the issue of back reaction due to the modified initial conditions. 



(25) 



(26) 



(27) 



B. Wigner function 

Another possible way of understanding the physical content of a quantum state, especially its classicality, is through 
the Wigner function (see, for e.g., Ref. [31]). Given a wave function V'k((Zk,'7), the Wigner function Wk (?k,Pk,??) is 
defined as [31] 



oo 

yVk (5k, Pk,??) = ^ j du] 



kV'k 



(«k + 2"k),»? 



V'k 



(Sk - 2"k),»? 



(28) 



The Wigner function corresponding to the Gaussian wave function (9) can be expressed as 



Wfe {qk,Pk,v) = - exp- 

TT 



iL + '^^(^)(fk-Jfe(??)9k)' 



(29) 
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where ak and Jk are given by 

al = {Rk + Rtr^ and Jk = i {Rk - R*k) ■ (30) 

On using the relations (A3) and (A5), we find that the quantities ak and Jk can be written in terms of the function 
fik and the Wronskian W{k) as follows: 

The quantum versus classical nature of the wave function depends on the evolutionary behaviour of (Jk{r]) and J^kiv)- 
It is possible for evolution to lead to a Wigner function sharply peaked around some region in the phase space, starting 
from a Wigner function which is uncorrclatcd in phase space [32]. For example, in the case of exponential inflation 
and the Bunch-Davies vacuum for the initial state which we will discuss later on (see Eqs . (43) and (55)), we we will 
find that 

-I = (^) (1 + kW) ; Jk = (^) (1 + k^iT' • (32) 



In such a case, Jkiv) — > 0, ct^ ^ oo at early times (ry — oo) which corresponds to a state sharply peaked around 
the g-axis. At late times (77 ^ 0), however, we have Jk{v) 00 with a finite cr^ which corresponds to a state sharply 
peaked around the p-axis. In fact, whenever fcry — > (corresponding to super Hubble scales), the Wigner function 
gets peaked around a classical trajectory. This can be verified more explicitly by studying the classical solution for 
our problem. Classically, for the case of exponential inflation, we can write the general solution for as: 



gk = 2Re 



* \ ikri 



(33) 



where C{k) is an arbitrary complex number. Writing C{k) — L{k) e^^'^^\ this solution becomes 



1L H 

gk = -2i Hr} cos \kr} + Uk)] H — sin [fcr? + Uk)] . (34) 

k 

The conjugate momentum pk = (a^ 9k) corresponding to the above gk is then given by 

Pi, = -{2Lk/H7]) iiin[kri + l{k)]. (35) 
The trajectory of the system in the phase space is given by 

At late times (when 77 0) or at super Hubble scales, we have (gk/Pk) ^ indicating a trajectory along the vertical 
p axis, which is precisely what we get from the Wigner function. On the other hand, one cannot naively take the 
early time limit (when r} — > —00) with finite gkjPk in the trajectory in Eq.(36) since gk becomes imaginary. One 
can, however, take the limit of 77 ^ — oo,pk keeping ?7Pk constant; in this limit, we obviously get a trajectory 
along the horizontal Pk = axis, which matches with the analysis based on the Wigner function. Incidentally, notice 
that the Hamiltonian for our system has a kinetic energy term K oc /a^ oc {prf)'^ and a potential energy term 
U oc c?q^ oc [q/rff. At late times, potential energy dominates over kinetic energy leading to near classical behaviour 
peaked aroimd g = 0; on the other hand, at early times if we let rj —00 keeping [pri) fixed, the kinetic energy 
remains finite and dominates over the potential energy. This is the quantum regime. (We plan to explore these ideas 
in detail in a separate publication). 

III. POWER LAW AND EXPONENTIAL INFLATION 

A. Standard initial conditions 

Let us now apply the above formalism to the case of power law inflation. Power law inflation corresponds to the 
situation wherein the scale factor a{t) grows with t as 

a{t)=aotP, (37) 
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where p> 1. In terms of the conformal time r], this scale factor can be written as 

a{v) = {-nv)^^+'\ 

where 7 and H are given by 



7 = - 



2p-l 
p-1 



and 



(38) 



(39) 



Note that 7 < —2 with 7 = —2 corresponding to exponential inflation. Also, the quantity H denotes the characteristic 
energy scale associated with inflation and, in the case of exponential inflation, it exactly matches the Hubble scale. 
On writing /Xfe = (/fc/a) in the differential equation (11), we find that fk satisfies the following equation: 



a 



fk 



0. 



(40) 



In a Friedmann universe described by the scale factor (38), the general solution to this differential equation can be 
written as (see, for e.g., Ref. [33], p. 362) 



ikrj) 



(41) 



where H^^^ and H^'^^ are the Hankel functions of the first and the second kind, respectively. The ^-dependent constant 
B{k) is to be fixed by choosing suitable initial conditions for each of the modes. For the above solution, it can be 
easily shown that the Wronskian condition (A5) leads to the relation (13) between B{k) and W(k). 

Let us first briefiy review the standard theory (see, for e.g., Ref. [28]) in which the initial conditions are imposed on 
sub-Hubble scales, i.e. when the physical wavelengths {k/a)~^ of the modes are much smaller than the Hubble radius 
. A natural choice for the initial condition will be the one in which each of the oscillators gk is in its ground state 
at sub-Hubble scales. This condition implies that the wave function (9) has the following asymptotic form: 



lim ^/ik (gk,??) 



{k/aH)—*oo 



cxp 



2« + 



(42) 



which, in turn, requires that, as (k/aH) 00, we need to have Rk -> {ka'^/2) and A^fe ^ (fcaVTr)^^" e'^'^''/^. These 



conditions can be satisfied provided — ^ {e^''^/y/2k) as {k/aH) 00 [cf. Eqs. (A3), (A4) and (A6)]. This can be 
achieved by setting B{k) = in Eq. (41), so that we have 



(43) 



and this choice corresponds to what is known as the Bimch-Davies vacuum [27]. Note that, according to the Eq. (13), 
B{k) = impHes W{k) = 1. Therefore, on substituting the above fk in Eq. (12) and imposing the condition of Hubble 
exit, viz. that (fc/a) = {zH), we obtain the spectrum of perturbations to be (for a recent discussion, see, for e.g., 
Ref. [34]) 



k'V^{k)=Ciz) [^] [to 



k 



2(7+2) 



where C{z) is given by 



In the limit of 2; — > 0, this expression simplifies to [34] 

'2-2(7+1) 



C(0) = 



2-JT 



i(7+i)^r 



7+ 



-(27+1) 



(44) 



(45) 



(46) 
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In the case of exponential inflation, corresponding to 7 = — 2, the perturbation spectrum in Eq. (44) reduces to 

t»p.(t) = (l±f!)^ (47, 

which is a spectrum that is exactly scale invariant. (Note that, for exponential inflation, Ti. = H .) As we pointed out 
before, the numerical value of the amplitude depends on whether we evaluate the expression at z = l,z = (27r) or as 

Let us now consider a more general situation with B{k) ^ 0. It is convenient to write 

B{k) = B{k) exp [i b{k)] , (48) 
so that the resulting power spectrum can be expressed as [on assuming that Hubble exit occurs at (k/a) = (zH)] 

/2 \ / L \ 2(7+2) 



fc3p*(fc) =C(^) (^^^ (^A^ ' [l-B^k)] ' (l + B^k)+ 2 B{k)cos[b{k)+ 7:^-29] 
where 6 is the phase of the Hankel function Hi^^ at Hubble exit, given by the relation 

If we write cos [b{k) + — 26] = d{k), then, the power spectrum (49) reduces to 



(49) 



(50) 



where — 1 < d{k) < 1. 

At super-Hubble scales (i.e. as 2 ^ 0), the above power spectrum bears a simple relation to the mean number of 
particles (n^) as given by Eq. (24). We find that they are related as follows: 

k^V.ik)^['^) (in,) ii. + Dzf) (^^(:^)'''"^ (52) 

This relation can be easily obtained by using the expression (25) for £k and Eq. (26) which relates £k to (rifc). At 
super-Hubble scales (i.e. as {kri) 0), the general solution for fk as given by Eq. (41) reduces to 

/fe(7?)(X-^(fc7?)^^+'^- (53) 

On using this expression, it is straightforward to show that /x^ = (/fc/a)' is sub-dominant to /Zfc at super-Hubble 
scales. Then, from Eqs. (25) and (26), we have 

^-^(IiIt) <-> 

It is then evident from the definition (12) that the power spectrum will be proportional to the mean number of particles 

at super Hubble scales. On using the above result, one can easily arrive at the relation (52), by first imposing the 
condition for Hubble exit [viz. that {k/a) = {zH)] and then taking the limit 2; ^ 0. In the limit of exponential 
inflation (i.e. as 7 ^ —2), we find that [fc^ 7^$(fc)] — {H'^ {nt) z'^/t^^) at super-Hubble scales. 

In the next section, we shall show that any modified spectrum obtained from a high energy model can be constructed 
in the standard theory by simply choosing suitable forms for the functions B{k) and b{k). But, before we do that, 
we shall discuss an alternative procedure for imposing the initial conditions wherein the initial condition for different 
modes are imposed at different times. 



B. Specifying initial conditions when [k/a{rik)] — Lp 



In the last section, we chose the quantum state by imposing a condition as (k/aH) 00 and one may question 
whether we are sure of the short distance, sub-Planck scale physics well enough to make this choice. An alternative 
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procedure, discussed extensively in the literature [16], attempts to address this question by choosing to impose the 
initial conditions for the different oscillators at different times. Specifically, one chooses the initial condition for the 
oscillator gk at a time r/k such that [k/a{rik)] — Lp^, where Lp is the Planck length. The initial state of the oscillator 
uniquely determines the function B{k) and the corresponding power spectrum can then be constructed using Eq. (49). 
However, it should be emphasized here that, apart from the fact that it has to be consistent with the observations, 
there exists no restrictions on the initial state to be chosen at rjk- 

To illustrate the point, let us consider the case of exponential expansion for which a{r]) = {—Tlr])~^. Let us 
assume that the oscillator is in the ground state at a given time, say, rji, which then requires that, at this instant, 
Rk ^ - {ia'^/2) (p'k/fJ.k) = (fcaV2). For the case of exponential expansion, the general solution for fk as given by 
Eq. (41) can be expressed in terms of simple functions as follows (see, for instance, Ref. [33], pp. 437-438): 



fkiv) = 



kri 



l + e^fc" + B{k) ( 1 - 7^ ) e-*'^'' 



kt] 



On imposing the condition Rk — {ka? /i) at r]i, we can determine the function B{k) to be 

B{k) = {l + 2ikri,y 



-1 g2ifc77i 



(55) 



(56) 



Our original choice of the Bunch-Davies vacuum corresponds to assuming rji — oo so that B{k) vanishes for all k. 
The modified procedure will be to choose 77, differently for different k by imposing the initial condition at, say, 
[k/a{r]k)] = Lp^. For exponential inflation, this translates to {krjk) = —{'HLp)~^ = a constant. For rji = rjk, 

the B{k) above reduces to [16, 20] 



S(fc)-[l-(2z/0]-'e- 



-(2i/{) 



(57) 



which is a constant independent of k. Thus the fc-dependence of the power spectrum remains unchanged when the 
initial condition on the k-th oscillator is imposed at a time such that [k/a{r]k)] = constant. The amplitude of the 
power spectrum, of course, gets scaled by a A;— independent factor; this is of no observable consequence since we do 
not know how to obtain the amplitude from a first principle theory anyway. 

The power spectrum corresponding to the above B{k) when evaluated at Hubble exit, say, when (k/a) = {z H), is 
given by 



k^r^ik) = 



1 + z' 



f7f_ 
V 27r2 



1 



2 



(z2-l)^2_4^^ 
2(l + z2) 



2z 



1 + Z2 



2z 



(58) 

This expression is exact and shows that the spectrum is strictly scale invariant, but is modified from the original result 
[viz. Eq. (47)] by a multiplicative factor independent of k. This factor depends on two parameters: (i) ^ = {H Lp) 
which measures the energy scale of inflation relative to the Planck scale and (ii) z which is the ratio of the Hubble 
radius and the physical wavelength of the perturbation. As we explained before, we can evaluate the power spectrum 
either at 2; = 1 or at 2 = (27r). For z = 1, the spectrum (58) reduces to 



k^ V^{k) 



while, for z = (27r), we get 



2f_ 

27r2 



2 



1 + 



^ cos 



e 



■ sm 



(59) 



(4^2,l)g2-87rg 
2 (1 + 47r2) 



cos 



(47r2 - 1) ^ + 2ttS; 
l + 47r2 



sm 



(60) 

The numerical value of these modified amplitudes will depend on the parameter ^ which is expected to be very small 
compared to unity for GUT scale inflation. For ^ <C 1, one can easily obtain the leading order terms of polynomial 
expressions involving ^, but determining cos(l/^) and sin(l/^) for ^ ^ 1 requires care. Since, one can easily replace 
(1/^) by, say, [(1/^) + (7r/2)], in the arguments of trigonometric functions to the leading order, these expressions are 
intrinsically ambiguous. For ^ <C 1 and z = 1, we get 



fc3 V^ik) ~ 



27r2 ; 



n- 



[l-^F{2/0], 



(61) 



where F is a rapidly oscillating cosine function. (Note that, if one evaluates the spectrum at super Hubble scales, i.e. 
as 2 ^ 0, then, instead of the cosine, one gets a sine function [16].) Similarly, when ^ <C 1, for 2 = (27r), we have 

■l + 47r2\ /7i:2\ / 47r^ \ . /(47r2 



27r2 



1 



47r^ 
l + 47r2 



l + 47r2 



G{2I0 



(62) 
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where G is a sine function and F, as above, is a cosine function. In each of these cases, a different choice for the 
sub-leading phase in the argument of the trigonometric function can make cosine into sine and vice versa. Hence, 
only the profile of the oscillating functions are of relevance in the limit of ^ <C 1, though the full expressions are often 
quoted in literature. 

Note that the quantum state with the choice of B{k) in Eq. (57) is a state with the initial condition for all the 
oscillators specified at a given moment of time. Thus the above analysis maps the prescription of specifying the 
quantum state for different oscillators at different times to specifying the initial condition at a given time. We can 
now explore the physical content of this quantum state at any given r] in terms of, for e.g., the mean occupation 
number (n^) in the instantaneous harmonic oscillator states. On using the expressions (21), (24), (55) and (57), we 
find the particle content of this quantum state to be 



{nk} 



1 



1 + y + fc^r?^ - I (^ - 4 fcr?) cos Q + 2fcr?j - ^ (1 + ^ fcr?) sin Q + 2kr] 



(63) 



This expression has several interesting features. 

First, let us consider very early times (i.e. as — > — oo) or very short wavelengths (i.e. as fc — > oo). It is precisely 
this limit which was considered uncertain due to trans-Planckian effects which motivated imposing the initial condition 
for different modes at different times; therefore, it is this limit which is of some interest to explore, to understand 
what kind of effective quantum state at ?7 = constant ^ — oo will lead to the vacuum state for the A;-th mode when 
(fcry) = — We see from Eq. (63) that the mean occupation number has two sets of contributions. The first three 
terms in Eq. (63) is secular, while the last two terms are oscillatory. The secular term increases monotonically from 
(^fe) = at ry = — oo. Thus, in the trans-Planckian limit (i.e. as 77 — > —00 or as fc ^ 00) all the modes have the 

same mean excitation (^^/4). (The oscillatory terms do not contribute in this limit.) 

Second, we do know that (rifc) must vanish for {kij) = since this is the condition we used to choose this state. 

As can be directly verified, this is indeed true for the expression in Eq. (63), but occurs because of a cancellation 
between the secular and the oscillatory terms. At later times, the secular terms dominate the oscillatory terms and 
(rife) increases on the average with superimposed oscillations. (We stress the fact that (rife) is computed in terms 
of the instantaneous harmonic oscillator modes; one should not think of them as particles which are produced with 
respect to the in-out states.) Thus, for each of the modes, we start with (n/j) = (^^/4) as ry — > —00, evolve to {uk) = 
at {kr]) = and grow to 



1 + + - I (^ + 4 0) cos f ^ - 2^) - e (1 - C ^) sin - 2z 



(64) 



at Hubble exit. We see that the alternative prescription of specifying the initial condition when the physical wavelength 
of the mode is comparable to the Planck length (which sounds reasonable at first sight) is equivalent to assuming a 
quantum state at a 77= constant hypersurface, with very specific properties. It can be realized only if unknown physical 
effects of the trans-Planckian sector acts in a particular manner to populate the modes with a specific prescription. 
It is far from clear whether this will be possible in a generic context. 

IV. MIMICKING THE TRANS-PLANCKIAN EFFECTS 

We shall now turn to the question of choosing an initial state such that a given power spectrum of perturbations is 
reproduced. Let us assume that, in the case of power-law inflation, a 'trans-Planckian' theory leads to the following 
spectrum: 

/ -xy2 \ / L \ 2(7-^2) 

[k^r^ik)]^ = Ciz)[^^j i^-j M{k), (65) 

where M(k) is the modification factor with, of course, M(k) > for all k. (Alternatively, let us suppose that a future 
CMB observation leads to such a spectrum, starting a hectic flurry of theoretical activity to explain it!) Evidently, 
the modifled spectrum (65) can be constructed from the spectrum (51) in the standard theory, provided we can find 
a positive definite function B{k) that satisfies the condition 

[l + B^{k) + 2B{k)d{k)] = [l-B'^{k)] M{k). (66) 

This is a quadratic equation in B{k) for each value of k and, if we choose d{k) = ±1, then, we find that the non trivial 
and positive definite roots for B{k) can be expressed in terms of M{k) as follows: 

= (67) 
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Given a M{k), this expression then allows us to reproduce the modified spectrum in the standard theory. However, 
note that the B(k) wc have constructed above is not unique as it corresponds to a particular value of the phase b{k) 
in Eq. (48) such that d{k) = ±1. Clearly, a whole class of such functions can be constructed, with different choices 
for d{k), but, as we shall see, this particular choice is sufficient for our purpose. 

As we had mentioned in Section I, the Planck scale modifications to the standard inflationary perturbation spectrum 
have been obtained in the literature by considering various high energy models for the scalar field [7] [20]. Most of 
these modified spectra deviate from the standard scale invariant spectrum at the ultra-violet end. On the other hand, 
the lower power in the quadrupole and the octopole moments of the CMB as measured by WMAP [6] — if it survives 
further releases of WMAP data!- -requires a suppression of power in the infra-red end of the spectnmi [35, 36]. Such 
spectra have also been obtained in certain high energy models of the scalar field [19, 22]. It is interesting to determine 
how these different modified spectra can be constructed by choosing suitable initial conditions on the scalar field in 
the standard theory. 

In the following subsections, we shall explicitly construct the function B(k) for four modified power spectra that 
have either been proposed to fit the observational data or have been obtained from a high energy model of the scalar 
field. 



A. Modified spectrum I 



Recently, the following spectrum was obtained in a power law inflationary scenario using a Lorentz invariant modifled 
theory [19]: 



[ev.{k)]^ = cii) (gj) (A 



2(7+2) 



where kc denotes the high energy scale, C is given by 



C= 2C(1)(7 + 1) 



3(7+1) 



(68) 



(69) 



and it is assumed that 10~^ < {Ti/kc) < lO"*^. (However, it should be mentioned here that this spectrum is not valid 
for arbitrarily small values of k [19].) In fact, a similar spectrum have also been obtained in non-commutative models 
of inflation [15, 22]. These spectra exhibit a suppression of power at the large scales, a feature that could be relevant 
to the low quadrupole moment observed in the CMB [35]. It is straightforward to see that, for z = 1, the function 
B{k) corresponding to the modified spectrum (68) is given by 



B{k) 



-(7+2) 



- 1 



(70) 



For 7 < —2, we have B — > A;('''+^) — > as fc — > oo. This implies that, while, towards the infra-red end, the initial state 
is different from the Bunch-Davies vacuum, the Bunch-Davies vacuum structure is retained at the ultraviolet end. 



B. Modifled spectrum II 



Another modified primordial spectrum that has been proposed (see, for e.g., Refs. [36]) in order to account for the 
low quadrupole moment observed in the CMB, is the following: 



[k^V^{k)\^ = Ak^-'-'^ 



-(k/ki)" 



(71) 



where A and rig are the scalar amplitude and index of the standard spectrum. The pivot scale k\ and the constant a 

(which turns out to be a positive number of the order of unity) are chosen to fit the CMB data. If we now assume 
that amplitude A and the index rig of the above modified spectrum are the same as those in the case of power-law 
inflation in the standard theory [cf. Eq. (44)], then, we find that B{k) is given by 



Note that B{k) -> 1 as A; 
end is not modified. 



B{k) 

and Sfe ~ e-('=/'=i)° - 



2g(fc/fci)° 



1 



(72) 



for large k. Once again, the vacuum structure at the ultraviolet 
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C. Modified spectrum III 



Another modified spectrum whose effects on the CMB has also been analyzed is the following spectrum [21] 



[fc^ P*(/e)] 



M 



2 7r2 



1-e 



sm 



2fk 



(73) 



where ^ ~ 10 e ~ 10 and ^2 is a pivot scale. In order to match the leading term in the spectrum, let us assume 



that the inflating background undergoes exponential expansion and that the spectrum is evaluated at z 
then straightforward to construct B{k) for the above spectrum and it is given by 



B{k) 



^{k/k2V |sin[(2/0(fc/fc2)^]| 
2-UfcA2)"'sin[(2/0 {k/k2T 



1. It is 



(74) 



D. Modified spectrum IV 

A more general modification of spectrum one can envisage is a spectrum which has corrections at both the infra-red 

and the ultra-violet ends. If wc now assmnc that the standard spectrum is modified in the infra-red as in Eq. (71) 
and has the same correction at the ultra-violet end as in Eq. (73), then the complete spectrum will be given by 



[fc^ Viik)] ^ = A A;("=-i) ( 1 - e-('=/'=i)° - ^ 



sm 



(75) 



where A and rig are the scalar amplitude and index of the standard spectrum and k\ and k2 are pivot scales such that 
ki <^k2- Also, as in the earlier case, let a be a positive constant of the order of unity. The corresponding B{k) can 
be easily obtained to be 



B{k)-. 



-{k/ki)" 



+ ^{k/k2)-'sm[{2/0ik/k2T 



2 _ e-Wki)" _ ^ {k/k2)~' sin [(2/0(^/^2)' 



(76) 



These examples demonstrate the fact that the modification of the spectrum due to 'trans-Planckian considerations' 
is degenerate with the choice of the initial quantum state. Without further input, such as an assumption for the 
choice of initial quantum state, observations cannot distinguish between these two physical effects. 



V. MODIFIED SPECTRA AND BACK REACTION 



An issue that remains unresolved in obtaining the modified spectra is whether the conditions that lead to modifica- 
tions of the standard spectrum will also lead to a large back reaction on the inflating background. In particular, will 

the energy in the quantum field dominate the inflaton energy thereby, possibly, terminating inflation? The approach 
we have adopted here allows us to address this issue along the following manner. 

Since the modified spectra from a high energy theory can be obtained from the standard theory with a suitable 
choice of initial conditions on the modes, we are probably justified in using these modes to evaluate the energy density 
of the quantum scalar field (for an earlier discussion on this point, sec the first reference listed under Ref. [12]). In 
what follows, we shall show that the energy density of the quantum scalar field above the Bunch-Davies vacuum is 
finite only if \M{k) — 1| decays faster than k~^ for large k and M{k) goes as k^ with — 2 < A < 2 for small k. Also, 
we shall restrict our attention to the easily tractable case of exponential inflation. 

Recall that the energy density of the quantum scalar field is given by Eq. (27) with the quantity £k to be evaluated 
using Eq. (25). In the case of exponential inflation, as we had mentioned earlier, the general solution for fk can be 
expressed in terms of simple functions [cf. Eq. (55)]. On substituting the solution (55) in the expression (25), we find 
that the energy density per mode of the quantum field is given by 

£k = (j^) [l - B^{k)] " ( [1 + B^{k)] {1 + 2 k^-q^) - 2B{k) cos[b{k) - 2krf\ + AB{k) [k-q) sin[6(A;) - 2kri\\ . (77) 
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In the Bunch-Davies vacuum, which corresponds to B{k) = 0, this expression for £k reduces to 

^'" = 1 + 2^ ™ 

and it known that the corresponding energy density (when regularized by point-splitting) is given by [27, 37] 

This energy density is much smaller than the energy density in the classical field that drives inflation. The energy 
density of the scalar field above the Bunch-Davies vacuum is then given by Eq. (27) with £k replaced by £k, where 

£k = £k- ^fe° = [1 - ^^W] (-S^(^) (1 + 2 A;^??^) - B{k) cos[6(fc) - 2krj\ + 2B{k) [kr]) sm[b{k) - 2kn]^ . 

(80) 

In order to understand the behavior of the above £k for large k, we can ignore the terms containing the sine and 
the cosine functions as they will oscillate rapidly in this limit. On neglecting these terms and, on making use of the 
relation (67) between B{k) and M{k), we find that, for large k, the above expression for can be written in terms 
of M{k) as follows: 

For the energy density p corresponding to this £k to be finite in the ultraviolet limit, the integral of (A;^ £k) over k 
should converge at large k. It is easy to see that this condition cannot be satisfied if M{k) oc k" for any value of v. 
For all values of this expression varies as at large k and hence the integral is divergent. To obtain a finite 

result, we need M{k) ^ 1 at large k. If we now assume that M{k) oc (l ± k~^) for large k, it is then clear that the 
energy density corresponding to the above £k will converge only if (5 > 2. Thus, the deviations from the standard 
spectra should die down faster than kr"^ for large k. 

One can also investigate the infrared limit in a similar fashion. As A; ^ 0, we find that the leading divergence arises 
due to the first term in the expression (80) for £k so that, in this limit, we have 

M{k)k ) ■ ^^^^ 

If we now assume that M{k) oc fc^ as fc — > 0, then the finiteness of the energy density p corresponding to this £k 

requires that — 2 < A < 2. It is therefore possible to enhance or reduce power in the infrared limit within a range 
and still maintain finite energy density. Modifications of the form [M(k) — 1] oc A^ with e > are also allowed and 
there are no restrictions on e in this case. (We stress that the finiteness of p is a necessary condition for ignoring 
back reaction, but it is not sufficient. The latter will require comparing the energy density in the quantum field with 
the background energy density which is difficult to do without assuming a specific model.) Clearly, amongst the four 
modified spectra that we have considered in the last section, only the second spectrum (provided a < 2) will lead to 
a finite energy density for the case of exponential infiation. 

Our expression in Eq. (81) shows that a constant M{k) independent of k leads to a divergent energy density above 
the Bunch-Davies vacuum. In particular, the state obtained by giving initial conditions for each of the modes at 
[fc/a(?7fe)]=constant (leading to a B{k) that is independent of k), produces a divergent contribution to the energy 
density and hence is suspect as a valid quantum state. It is sometimes argued in the literature that this state has the 
same energy density as the Bunch-Davies vacuum along the following lines: When B is independent of fc, we are dealing 
with mode functions of the form ( fk + B f^), with fk given by Eq. (55), which belong the set of so called a- vacua [38]. 
It is possible to construct a regularization scheme such that the divergences which arise at the coincidence limit of 
the Greens function in this case, is the same as that in the case of ;B = 0. Such a subtraction will lead to p = pbd 
for these states. In our approach, this is equivalent to ignoring the contributions in Eq. (81) when M{k) = Mq is a 
constant different from unity. Since this leads to £k oc k, one can think of Mq as the (constant) occupation number 
(rifc) for all k; the regularization involves ignoring all these 'particles' in measuring the energy density. 

We believe this argument is spurious for two reasons. Firstly, a transparent and direct discussion in terms of 
harmonic oscillators presented above gives a different result showing that M{k) = Mo= constant leads to a divergent 
energy density. So, at the least, the results depend on the explicit regularization procedure adopted. Secondly, in 
a realistic model, we will have to deal with a weak dependence of B on rather than no dependence. Then, the 
argument based on a-vacua will no longer hold, irrespective of how weak this dependence is. The analysis given 
above, however, is completely transparent and clear. 
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VI. DISCUSSION 



The fact that the measurements of the CMB anisotropies strongly indicate a primordial spectrum that is nearly 
scale invariant already implies that the initial state of the quantum scalar field is a state that is 'very close' to the 
Bunch-Davies vacuum. Clearly, sufficiently precise measurements of the anisotropies in the CMB can provide us with 
the form the inflationary perturbation spectrum to a good accuracy. However, this information can at most help 
us determine the initial state of the quantum scalar field in the standard theory and it is not sufficient to aid us in 
discriminating between the various Planck scale models of matter fields. The reason being that, given a perturbation 
spectrum, we should be able obtain the spectrum from any high energy model of the matter field by simply choosing 
a suitable state in the modified theory just as we had done in the standard theory. 

The above result can be obtained either in the Heisenberg picture or in the Schrodinger picture since these de- 
scriptions are mathematically equivalent. We have, however, adapted the Schrodinger picture since the problem we 
are discussing is essentially that of a harmonic oscillator with a time dependent frequency for which the intuition 
available in the Schrodinger picture is of some value. Our discussion clearly shows that virtually any power spectrum 
which is either observed in future or suggested by phenomenological models, can be reproduced by a suitable choice of 
the quantum state. Evidently, without additional assumptions one cannot disentangle the dynamics from the initial 
conditions and, hence, the CMB observations alone cannot act as a discriminator between different theoretical models. 

There are several possible avemics for future work arising from this discussion. One particularly interesting question 
will be the evolution of the quantum state of the universe into the future. Several recent observations [39] (as well as 
not so recent observations, see Refs. [40]) suggest that the universe has just entered an accelerating phase dominated 
by dark energy with an equation of state P « —p. While the nature of this dark energy is unclear, it seems likely that 
at least at sufficiently large scales it will act like a cosmological constant leading to a late time de Sitter phase [41]. 
It will be interesting to study the late time evolution of the quantum wave function of the scalar field. The de 
Sitter phase in the future should lead to its own thermal fluctuations with a characteristic temperature and it will be 
interesting to see how that can emerge. 

Another issue of interest is the study of correlations across the horizon in the case of de Sitter spacetime. It has 
been shown that the quantum entanglement of modes across the horizon can lead to a holographic description of 
gravity [42] . The effects of this entanglement are easy to calculate when de Sitter spacetime is described in the static 
coordinates. On the other hand, the time dependent Gaussian state used in this paper is more naturally tuned to 
the Friedmann coordinates of the de Sitter spacetime. It will be worthwhile to relate these two descriptions and 
understand how the correlations across the horizon arises in the time dependent background. 

Finally, the description in terms of the wave function can be easily translated to one in terms of the path integral 
kernel using the Feynman-Kac formula. This will allow one to provide a purely path integral derivation of the results 
presented in this paper. All these issues are currently under investigation. 



APPENDIX A: EVALUATING THE POWER SPECTRUM 



Recall that the quantum state Vk of the time dependent oscillator was described by the Gaussian wave func- 
tion (9). On substituting the wave function (9) in the Schrodinger equation (7) and equating the coefiicients of qu 
and q^, we find that Nk and Rk satisfy the following differential equations: 

iiv; = (Ai) 



where the prime, as before, denotes differentiation with respect to r]. Let us now introduce a new quantity /ife(r?) 
through the relation [28-30] 

On substituting this expression in the above differential equation for Rk, we find that fj,k satisfies the differential 
equation (11) which is the same as the classical equation of motion satisfied by the classical oscillator variable gk- 
Also, in terms of Hk, the differential equation (Al) can be integrated to obtain 

*^(^), (A4) 
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where D{k) is a fc-dependent constant determined by the normahzation condition (10). Note that the differential 
equation satisfied by ii^ [viz. Eq. (11)] imphcs the Wronskian condition 

il^kl^'^* -l^'kl^l)=-[iWik)/a''], (A5) 

where W{k) is a fc-dcpendcnt constant. The relations (A3) and (A4) along with the above Wronskian condition and 
the normalization condition (10) determine D{k) to be 

Thus, the solution to classical equation of motion allows us to construct the wave function for the corresponding 
quantum problem [25, 28-30]. However, it should be pointed out here that, while is real, /Xfc, in general, is a 
complex quantity. 

The power spectrum of the perturbations (8) is given by 

oo oo 

ev^^k) = ^ J dqu\^^\'ql='^^^ J dquqle-i^^^+^'^^^i], (A7) 

— CO — oo 

where we have substituted the expression (9) for the wave function ip^. On carrying out the integral over in the 
above equation and making use of the relation (10) and the Wronskian condition (A5), we obtain the power spectrum 
to be 

which is expression (12) we have quoted in the text. 
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